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Abstract 

We provide a problem definition of the stable marriage problem for a general number 
of parties p under a natural preference scheme in which each person has simple lists 
for the other parties. We extend the notion of stability in a natural way and present 
so called elemental and compound algorithms to generate matchings for a problem 
instance. We demonstrate the stability of matchings generated by both algorithms, as 
well as show that the former runs in 0{pn‘^) time. 


1 Introduction 

The stable marriage problem (SM) is a famous problem in mathematics in which 
there exists a community of n men and n women, all of whom are to be paired with each 
other heterogeneously in marriage. Each individual provides a complete preference list 
ranking the members of the opposite sex according to his or her preference for marriage. 
The final matching sought is one in which no two people would rather be married to 
each other over their current spouses, and is thus called stable. It is of note that 
stability is a heuristic approach to optimizing matching. 

In 1962, David Gale and Lloyd Shapley [2] presented an algorithm to solve SM. 
In the Gale-Shapley (GS) algorithm, each man proposes to his favorite woman, and 
each woman is temporarily matched to the man that proposes to her. If a woman 
is proposed to by multiple men, she is temporarily matched to the proposing man 
she prefers the most. If there remain any unmatched people after the first round of 
matching, the single men propose again, this time to their second choices. If a woman 
who has already been matched is proposed to by a man whom she prefers over her 
current partner, she will leave her current partner and become paired with the new 
man. This process is repeated, with each unmatched man proposing to the women in 
order of his preference list and each woman choosing her best possible mate who has 
proposed to her, until everyone has been matched. 

In [2], the following properties of the GS algorithm were shown: 

• It terminates 

• It is O(n^); maximal number of rounds is — 2n + 2 

• Resultant matching is stable 

• matching is optimal (of stable) for proposing party 
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• matching is pessimal (of stable) for responding party 

For a thorough introduction to the stable marriage problem, including the Gale- 
Shapley algorithm, we refer the reader to [3]. 

In 1976, Donald Knuth proposed twelve open questions j6j on SM, one of which 
asked to generalize SM from two to three parties—the 3-dimensional SM (3DSM). 
Given the open-ended wording, the question has been addressed under various inter¬ 
pretations with respect to structure of preferences and dehnitions of stability. With the 
addition of indifference in preferences [5] , there have been many different constructions 
of 3DSM, most of which have been shown to be AAP-complete [Z][l]. 

However, 3DSM has been shown to work under a simple scheme [T] in which each 
individual provides two simple preference lists for the other two parties. Along this 
vein, this paper will consider the p-dimensional SM (pDSM), where p>2, and propose 
two types of algorithms to deal with it. Before presenting the algorithms in their 
entirety, we motivate them with the simplest nontrivial consideration, when p = 3. 

In an instance of 3DSM, we have a community of men, women, and dogs. One 
possible way to match the men, women, and dogs together is to choose two parties, 
say men and women, and create a matching between them using the Gale-Shapley 
algorithm. 

Repeating this with another pair of parties, say women and dogs, we arrive at sets 
of man-woman and woman-dog pairs, from which we can deduce a matching for the 
entire community. This is the idea motivating the elemental algorithm. 

Alternatively, we can take the man-woman pairs and view each as members of a 
compound party, the humans. Each human’s preferences are constructed by combining 
each (associated) man’s and woman’s preferences for dogs. Similarly, we can modify 
each dog’s preferences by combining its preference of men and women into preferences 
for humans. This is the idea motivating the compound algorithm. 

Here the names elemental and compound are chosen in allusion to chemistry; el¬ 
ements are the fundamental building blocks, and compounds are created by bonding 
elements together. In the elemental algorithm, each party is treated as an individual 
element in “pure” form, whereas in the compound algorithm, parties are “bonded” 
together via stable matchings. 


2 Problem Definition 

We proceed with a formal construction of the multidimensional stable marriage 
problem. An instance of the p-dimensional stable marriage problem (pDSM) is an 
ordered pair {V, L) with a set V of p disjoint parties with n elements each, and a set 
L of preference lists for every element (to be defined below). An element of a given 
party is said to be a member of that party. Let Up = IJpg-p P be the community 
of V. For all x € Up, let prt{x) return the party of which x is a member. 

Associated with each x G Up is a p — 1 x n strictly ordered preference array 
Lx of x’s preferences defined as follows. For all y G P € P where x ^ P, let Lx{y) = j 
denote that y is x’s preferred member of P. Let 

Lx{P) = {y G P I strictly ordered by Lx{y)} 

Lx = {Lx{P) I P G P} 
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L — {Lx I X € Uj)} 

Furthermore, given a,b € P, let a >-x h denote that Lx{a) < Lx{b) and let a Px b 
denote that Lx{a) < Lx{b). 

A family F C U-p is a set of p elements, one member from each party. A matching 
P over P is a partition of Up into n families. Elements of a single family are said to 
be relatives in P. Let relp{x, P) return x’s relative in P from P. 

For a given P, a family F ^ P is blocking if and only if 

• X Py relp{y,prt{x){ , y Px relp{x,prt{y){ for all x,y € F 

• for each x, there exists z & F such that z Px relp[x,prt{z)) 

A matching P is unstable if there exists a blocking family in Up. P is otherwise 

stable. 

There are several possible ways to define the problem, with regards to both manner 
of describing preferences and definition of stability. The manner of preferences was 
chosen, in part, because of its simplicity. In 3DSM, simple preference lists leads to an 
efficient algorithm, while more permutation-oriented (Cartesian product) setups have 
generated complex problems proven to be A/'P-complete. The definition of stability 
chosen resembles that of the traditional problem most closely. This is because if there 
is a blocking family of elements who all prefer each other to their corresponding partners 
in T”, they would all “elope” and desert their established families in real life. Whereas 
setups that permit indifference among individuals have led to A/'P-complete problems 
as well [5]. 

Two types of algorithms will be presented for pDSM, elemental and compound, 
both of which are novel extensions of the original GS algorithm. Given any algorithm 
A, let (A) denote the matching generated by A. Given a set of algorithms A, let 
(A) = {(A) I A G A}. 

3 Elemental Algorithms 

A tentative definition of an elemental algorithm will initially be provided, which will 
prompt a deeper understanding and a more rigorous definition. Let GS{P,Q) denote 
that P proposes to Q according to the GS algorithm. An important way of viewing 
GS{P, Q) is that it establishes a bijection between P and Q. 

An elemental algorithm is a set £ of bijections GS{R, S) such that for all P,Q 
there exists a unique bijection (either directly, or indirectly by composition) between 
P and Q. To execute e, each element in e is executed, generating a unique 1-1 corre¬ 
spondence between each pair of parties, and thus a matching for the problem. 

In SM, the GS algorithm is executed on the pair of men and women. This may 
be viewed as a directed graph with a vertex for each party and an edge directed from 
the proposing into the responding vertex. Unless otherwise specified, all graphs are 
assumed to be simple and labeled. 

In pDSM, these graphs may be similarly constructed. Given an elemental algorithm 
£, an elemental graph G = {V,E} is generated according to the following bijections. 

V^V :P^Vp e^E: GS{P, Q) ^ e^q 

The following theorem will provide a basis for the formal definition of an elemental 
algorithm. 
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Theorem 1. G is an elemental graph G is a tree. 

Proof. This bidirectional proof will be split into two parts. The first will be proven by 
contradiction, and the second by direct construction. 

part 1: Take an elemental graph G. Since Cpq represents a bijection between P and 
Q, an (undirected) walk W 'm. G establishes a composition of bijections, relating all 
parties with corresponding vertices in W to each other. Since a bijection is established 
between all pairs of parties, G must be connected. Further, since each bijection is 
unique, there can be no cycles in G. So G is a tree. 

part 2: Take a tree T = {!/, FI}. Let T generate a set 5 of bijections GS{P,Q) 
according to the following bijections: 

V ^V:Vp^P E ^8: epqe^ GS{P, Q) 

Since T is a tree, there exists a unique walk between Vp and Vq, corresponding to a 
unique bijection between P and Q. Thus, by definition, 5 is an elemental algorithm. 
Since the maps used to generate 5 from T were bijective, 8 may be thought to generate 
T using the inverse maps. Then, since 8 is an elemental algorithm, T is an elemental 
graph. □ 

The formal definition follows from the reasoning in part 2 of Theorem 1. Given a 
directed tree T = {V, E} with p vertices, an elemental algorithm st over P is a set 
of bijections GS{P,Q) generated according to the following bijections: 

V ^ V '. Vp ^ P E ^ tT '■ Vpq 1-^ GS{P, Q) 

The execution of st remains unchanged; each element of eT{E) is executed, generating 
a matching 

It is important to note here that a directed tree is specified. The direction of the 
edge indicates which party proposes, an important fact when executing the algorithm. 
However, such specification was unnecessary in the proof of Theorem 1 because both 
are viable elemental algorithms. 

Having provided a definition, we now show that the elemental algorithm shares 
many of the nice properties the GS algorithm has. Let Tp be the set of all undirected 
trees with p vertices. Let S{P) = {exiE) \ T G Tp}. For ease (given appropriate 
context), let e, T and £ be short for zt, Tp, and £{V), respectively. 

Theorem 2. All elemental algorithms terminate. 

Proof. Take e G F. By construction, z contains p — 1 instances of the GS algorithm. 
Since the GS algorithm terminates [2] and p — 1 is finite, e must terminate. □ 

Theorem 3. AH elemental algorithms yield stable matehings. 

Proof. Consider an arbitrary e G F. Suppose (e) is unstable; there exists a blocking 
family E ^ (e). Take x,y ^ E where x G P,y G Q. 

case 1: GS{P,Q) G z. Since GS{P,Q) is stable [2], we have x = rel(^^'^{y, P). 
case 2: GS{P,Q) ^ z. Take the tree T = {V,E} that generated z. There is a 
unique walk VF in T from Vp to Vq. By relabeling vertices, without loss of generality, 
let 

W = U 1 U 2 U 3 ■ ■ - Uw for some 1 < w < p, where ui = Vp and Uw = Vq 
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Relabeling parties accordingly, we have that 

GS{Pk, Pk+i) G £ for 1 < k < w, where Pi = P and Pw = Q 

We now apply case 1 w — 1 times on Xfc,Xfc+i, where Xk € Pfc,xi = x and Xy^ = y. 
This gives Xk = (x^+i, Pfc). Then by transitivity, xi = rel(^^'f{xw, Pi), or x = 

rel{^){y,P). 

Combining both cases, we have that P G (e), which is a contradiction. Therefore 
(e) is stable. □ 

For 3DSM, this setup was criticized for its simplicity: “It is not hard to see that 
we can apply the Gale-Shapley algorithm twice to get a weak stable matching: letting 
the men propose to women and then propose to dogs. Women and dogs make the 
decision of acceptance or rejection based on their simple lists of men.” [3] However, 
previous literature has overlooked the multitude of combinations of matchings available, 
allowing one to customize the algorithm to fit the given task at hand. The following 
result computes the exact number of such options. 

Theorem 4. \£\ = 2P~^pP~‘^ 

Proof. Since each T £ T generates a distinct e{T) G £, we have that |T| = \£\. By 
Cayley’s Formula, there are undirected trees with p vertices. Since there are p — 1 
edges to a tree, each undirected tree corresponds to 2^“^ distinct directed trees. Thus 
\£\ = \T\ = 2P-V-2. □ 

Theorem 5. An elemental algorithm is 0{pn^). 

Proof. Recall that the GS algorithm is O(n^) and the maximal number of rounds is 
— 2re + 2. Since an elemental algorithm applies the GS algorithm p—1 times, the 
maximal number of rounds in an elemental algorithm is (p — l)(n^ — 2n + 2). Therefore, 
it is 0 {pn^). □ 

4 Problem Structure 

Before dehning compound algorithms formally, we need to develop a vocabulary for 
some of the natural structure of a multidimensional stable marriage instance. 

Given a pDSM P, a (?DSM Q is said to be a subproblem of V if and only if Q C P. 
Given P, a set partition 7r(P) ^ P, {P} is said to be a problem partition containing 
subproblems of P. Let n(P) be the set of all problem partitions. For ease, let vr, H be 
short for vr(P), n(P), respectively. 

Given a problem partition 7r(P) containing 1 < p' < p subproblems with a matching 
P over each Q G vr, then a reduced problem (P', L') is a p'DSM where V = tt and 
L' is defined as follows. Given matchings P, Q for subproblems Q,TZ £V', respectively, 
for all families F £ P,G £ Q, let 

L'p{G)= Y. 

xGF, y£G 

Lp{TZ) = {G £ Q I strictly ordered according to Lp{G)} 

L'p = {Lf(P) I P G P'} 
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L' = {L'p I F G t/p/} 

The definition of the reduced problem effectively collapses each subproblem Q a V' 
into a single party and each family F £ F into a single individual. 

Given a matching F' over reduced problem {V\F), F' may be expanded to F 
according to 

/■= U F' 

F’&T’ 

to give the equivalent matching F over the original problem (F, L). 


5 Compound Algorithms 

Given a pDSM (F,L), a compound algorithm C is executed over V according to 
the following recursive procedure. Two counters i, c are both intially set to zero. 

1. Take an elemental algorithm e G £{F). If n(F) ^ 0, take a problem partition 
TT G n(F). Else, let (C) = (e) and go to step 6. 

2. Greate reduced problem set V' = vr; for each subproblem Q G vr, take 

an elemental algorithm e G S{Q) and generate L' using matchings (£(Q)}. 

3. Index i by one. 

4. If |7r| > 1, repeat from step 1 letting {V,L) = (V',L'). Else, let c = i. 

5. Given {z{V')) used in the execution of step 2, construct {C) as the result of 
expanding {ziV')) c times. 

6. Return {C). 

Having provided a definition, we now show that compound algorithms share many of 
the nice properties the elemental algorithms have. Let C('P) be the set of all compound 
algorithms over V. 

Theorem 6. All compound algorithms terminate. 

Proof. The statement will be proven by induction on |F|. Take C G C{V). 

Base case: \V\ = 2 

Let V = {P,Q}. Without loss of generality, £ = GS{P,Q). The only possible set 
partitions are V, {V}, neither of which are problem partitions. Thus n(F) = 0; for 
the execution of C, we go from step 1 to step 6, return {GS{P,Q)), and terminate. 

Inductive step: Given {V] = k + 1 and C{Q) terminates when \Q\ = i for all 
2 < i < k. 

Recall that a problem partition omits the trivial partition vr = {V}. This ensures, 
for all TT G n('P), that \V\ > |Q| for all Q G tt. Thus k + 1 > |Q|, meaning that \Q\ = i 
for some 2 < i < k. By the inductive hypothesis, G{Q) terminates. Therefore, so does 
G{V). □ 

Before showing stability, we need a lemma. 

Lemma 1. Given pDSM {V,L) reduced to {V',L') using stable matchings Fq for all 
Q £ V. If matching F' over V' is stable, then the expanded matching F over V is 
stable. 
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Proof. The proof will be given by taking any Q,TZ (zV and showing that the expansion 
of T”' n (Q U 7^), which is TTl (Q U 7^), is stable. For ease of notation within the proof 
itself, assume that F' stands for T”' fl (Q U 7^), and F for T” D (Q U 7^) . 

Take Q, 7^ € V'. Let Fq,Ftz be the stable matchings over Q, 7^, respectively. 
Suppose F is unstable; there exists a blocking family F ^ F. Since F is blocking 
and F' is stable, F cannot be an expansion of some F' G F'. Therefore, since Fq^Ftz 
are stable, F = (G n Q) U {H n TZ), for some G,H G F'. 

For ease of notation, let 

Fi = G n Q, Fa = G n 7^, F 3 = F n Q, F 4 = F n 7^ 

Thus, without loss of generality, F = Fi U F 4 . By definition of V', we have that 
Fi,F2,F3, F 4 G Fp/. 

Consider x G Fi,y G F 4 . Since x,y € F and F is blocking, we have that x Fy 
reljz[y,prt{x)) and y reljz[x,prt{y)). This means that 

Lxiy) < Lx{reljz{x,prt{y))^ 

Ly{x) < Ly[reljz[y,prt{x))^ 

Summing over all such x, y pairs, 

X] Lx{y)< Fx(relr{ x,prt{y)^'^ 

xeFi,yeF4 x£Fi,yeF4 

Lyix)< Y Ly{relr[y,prt{x))^ 
x&Fi,y£F4 x£Fi,yeF4 

Additionally, for each x G Fi,y G F 4 there exist z G F 4 ,w G Fi such that z Fx 
reljz{x,prt{z)) and w >-y reljz{y,prt{w)). This implies that 

X] Lx{y)< Y x,prt{y))'^ 

x&Fi,yeF4 xeFi,yeF4 

Y < Y Ly(^reljz{y,prt{x))'^ 

xeFi,yeF4 x£Fi,yeF4 

Then, by definition of L', we have 

L'p^(F4) < L'p^(F 2 ) 

L'p^(Fi) < L'p^(F3) 

But this contradicts the fact that F' is stable. Therefore there can be no blocking 
family in F, and thus F is stable. □ 

Theorem 7. All compound algorithms are stable. 

Proof. This statement will be proven by tracing the steps of the compound algorithm 
and proving that stability is preserved throughout. Take G G C{V). Upon executing 
the algorithm, let c be the final number of times that steps 1 — 4 were repeated. 

Nothing occurs with regards to stability in steps 1, 3, 4, and 6 . In step 2, by 
Theorem 3 each matching (e(Q)) is stable for all Q G vr. In step 5, (e(F')) is stable by 
Theorem 3. The initial algorithm conditions and steps 2 and 5 satisfy the givens for 
Lemma 1, so by applying Lemma 1 c times on {e(V')), we have that the expanded 
matching of (e(F')), (G), is stable over the original V. □ 
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6 Open Problems 

There are several interesting questions we can now ask. 

• determine optimal and pessimal matchings for a given party 

• determine egalitarian matchings 

• determine efficiency of the compound algorithm 

• determine additional structure of matchings under elemental and compound al¬ 
gorithms (dependent on chosen directed tree, etc.) 

• find applications 
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